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Estimation of the mass outflow rates from viscous 
accretion discs 
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ABSTRACT 

We study viscous accretion disc around black holes, and all possible accretion solu- 
tions, including shocked as well as shock free accretion branches. Shock driven bipolar 
outflows from a viscous accretion disc around a black hole has been investigated. One 
can identify two critical viscosity parameters a c i and a cu , within which the stationary 
shocks may occur, for each set of boundary conditions. Adiabatic shock has been found 
for upto viscosity parameter a — 0.3, while in presence of dissipation and massloss 
we have found stationary shock upto a = 0.15. The mass outflow rate may increase 
or decrease with the change in disc parameters, and is usually around few to 10 % of 
the mass inflow rate. We show that for the same outer boundary condition, the shock 
front decreases to a smaller distance with the increase of a. We also show that the 
increase in dissipation reduces the thermal driving in the post-shock disc, and hence 
the mass outflow rate decreases upto a few %. 
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1 INTRODUCTION 

One of the curious thing about the galactic as well as ex- 
tra galactic black hole candidates is that, they show mod- 
erate to strong jets. Since black holes are compact and do 
not have hard surfaces, outflows and jets can originate only 
from the accreting material. And observations have indeed 
showed that the persistent jet or outflow activities are cor- 
related with the hard spectral states of the disc i.e., when 
the power radiat ed maximizes in th e hard non-thermal part 
of the spectrum (|Gallo et. alj|2003h . Moreover, there seems 
to be enough evidence that these jets are pro duced within 
100 S chwarzschild radius of the central object jjunor et all 
1 19991 ). So accretion disc models which are to be considered 
for persistent jet generation should allow the formation of 
jet close to the central object, and that the spectra of the 
disc should be in the hard state. 

One of the most widely used accretion disc model 
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disc was chopped off adhoc, and cannot automatically 
explain jet generation mechanism, neither can it explain 
the relative small size of the jet base. Hence, one has to 
consider accretion models which have significant advection. 
An elegant summary of all typ es of viscous acc retion disc 
model has been provided by iLee et all (|201li ). For the 
sake of completeness, we present a very brief account of 
advective discs, in the following. 

Accretion onto black hole is necessarily transonic since 
radial velocity can only be small far away from the black 
hole, but matter must cross the horizon with velocity equal 
to the speed of light (c). Hence accretion solutions around 
black holes should consider the advection term self con- 
sistently. Accretion solution like Bondi fl ow or, radial flow 
|Bondll952l ; IChattopadhvav fc Rvul2009l ) although satisfies 



the Kepleria n disc 



(jShakura 

iNovikov fc Thornd 1 1 9731 1. is very successful in explain- 
ing the multicoloured black body part of the spectrum. 
However, the presence of Keplerian disc alone is inadequate 
in explaining the existence of the non-thermal tail of the 
spectrum. Moreover, the advection term in the equation of 
motion is poorly handled in this model, where the inner 
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the transonicity criteria, but is of very low luminosity. There- 
fore, there was a need to consider rotating flow which are 
transonic too, such that the infall timescale would be long 
enough to generate the photons observed from microquasars 
and AGNs. The accretion m odel with adve c tion which got 
wide attention, was ADAF (|lchimarulll977l ; iNaravan fc Yil 
Il994f ). Initially ADAF was constructed around a Newto- 
nian gravitational potential, where the viscously dissipated 
energy is advected along with the mass, and the momen- 
tum of the flow. The original ADAF solution was self- 
similar and wholly subsonic, thus violating the inner bound- 
ary condition around a black hole. This inadequacy was re- 
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stored when global solutions of ADAF in presence of strong 
gravity showed that the flow actually becomes transonic 
at around few Schwarzschild radii (r 9 ), while remains sub- 
sonic elsewhere. The self-similarity of such a s olution may 
be m aintained far away from the sonic point |Chen et all 
1 19971 ). Moreover, the Bernoulli parameter is positive mak- 
ing the solutions unbound. This prompted the introduc- 
tion of self-similar outflows in the ADAF abbreviated as 
ADI OS, which would carry away m ass, momentum and en- 
ergy (Blan dford fc Begelmanlll999r ). Although this is an in- 
teresting addition to the generalization of ADAF type so- 
lutions, no physical mechanism was identified except the 
positivity of Bernoulli parameter of the accretion flow, that 
would drive these outflows. 

Simultaneous to the research on ADAF class of so- 
lutions, a lot of progress was being made in the re- 
search of general adv ective, rotating sol u tions. For ro- 
tating advective flow, iLiang fc Thompsonl (|l980T l showed 
that with the increase in angular momentum of the invis- 
cid flow, the number of physical critical points increases 
from one to two, and later it was shown that a stand- 
ing shock can exist in between the two critical points 
(lFukuell987l;IChakrabartill 19891; iFukumura fc Kazanasll2007i ; 
IChattopadhvavll2008l ; IChattopadhvav fc Chakrabartill201lT ). 
Since the accretion shock is centrifugal pressure mediated, 
there was apprehension about the stability and formation 
of such shocks, in presence of processes such as viscosity, 
which transports angular momentum outwards. All doubts 
about the stability of su ch shocks in presence of viscosity was 
subsequently removed (|Lanzafame et q/1ll998l; IChakrabartil 



IChakrabarti fc Titarchukl dl995h; IChakrabarti fc Mandall 
|2006MMandal fc Chakrabartil |2008l . l2010l ) showed that the 
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Das fc Chattopadhvayll2008l; ILee et 0^.112011^ 7 Furthermore. 
Fukumura fc Tsurutal ((200J) conjectured about the exis- 
tence of multiple shocks and later, independent numerical 
simulations serendipitously found the existence of transient 
multipl e shocks in presence o f Shakura Sunya ev type vis- 
cosity jLanzafame et aZ.I I2OO8I ; iLee et aEll201ll ). Although, 
both general advective and ADAF solutions start with the 
same set of equations, it was intriguing that there can be 
two mutually exclusive class of solutions, and without any 
knowledge und er which condition these solutions may arise. 
iLu et all\ 19991 ) later showed that the global ADAF solution 
is a subset of general advective solutions. In other words, 
the models which concentrate only on the regime where 
the gravitational energy is converted mainly to the ther- 
mal and rotational energy, may either be cooling dominated 
(e.<7.,Keplerian disc; Shakura & Sunyaev, 1973; Novikov 
fc T horne, 1973) or advection dominated (|Naravan et all 
1997). Either way, these models remain mainly subsonic 
(except very close to the black hole) and hence do not 
show shock transition. However, if the entire parameter 
space is searched, one can retrieve solutions which have 
enough kinetic energy to become transonic at distances of 
few x 100r 9 . A subset amongst these solutions admit shock 
transitions when shock conditions are satisfied. This is the 
physical reason why som e disk models show shock trans i- 
tions and others do not |Lu et al. I ll999l ; iDas et al.ll2009l ). 
Whether a flow will follow an ADAF solution or some kind 
of hybrid solution with or without shock will depend on the 
outer boundary condition and the physical processes domi- 
nant in the disc. 

The shock model was later used to explain observations. 



post-shock region being hotter, can produce the hard power- 
law tail by inverse-Comptonizing soft photons from pre- 
shock and post-shock parts of the accretion disc. The soft 
state and hard states are automatically explained depend- 
ing on the existence or non-existence of the shock. Pres- 
ence of such transonic advective flow has also been sug- 
gested by observations |Smith et a/.ll200ll . I2002T ). Infact the 
'hardness- intensity-diagram', or HID a hysteresis like be- 
haviour se^riJri_jniCTO^iiasarehas_been reproduced by this 
model (|Mandal fc Chakrabartil l20ld ). In other words, the 
post-shock region is the elusive Compton cloud. 

Interestingly enough, the shock fronts are stable in 
a limited region of energy-angular momentum parameter 
space and naturally give rise to time dependent solutions 
whenever the exact momentum balance across the shock 
front is not achieved. This mi ght be due to different rates of 
cooling dMolteni et a/.lll996b|). or different rates of viscous 
transport (jLanzafame et aZ.lll99Sl : ILee et aZ.ll201ll ). These 
oscillations we re also confirmed by pure general relativis- 
tic simulations jAoki et a/.l2004l ; lNagakura" fc Yamadal2008l . 



120091 ). iMolteni et all (|l996bfi ~~suggested that if the post 
shock region oscillates, then the hard radiation produced 
by the post-shock region would oscillate as well, and hence 
explain the QPO. Infact, the evolution of QPO frequen- 
cies during the outburst states of various microquasars like 
XTEJ1 550-654, GRO 1655-40 etc were explained by this 
model l|Chakrabarti et 0/1120091 ) by assuming inward drift of 
the oscillating shock due to increased viscosity of the flow. 

Another interesting consequence of accretion shock is 
that, it automatically explains the formation of outflows. 
The unbalanced pressure gradient force in the axial direction 
drives matter in the form of outflows and may be considered 
as the precursor of jets dMolteni et all 1994 1_ Chakrabart 



1999 1; IChattopadhvav fc Das! 120071 ; IDas fc Chattopadhyay 
200a ) . These outflows can be accelerated by various 



accelerating mechanisms t o relativistic terminal speeds 
IChattopadhvav et all |2004| ; IChattopadhvav! 120051 ). If the 
shock conditions are properly considered, then the mass 
outflow should reduce the pressure of the post-shock re- 
gion and hence the shock would move inward, which in 
turn, would modify the shock paramete r space in presence 
of massloss jSingh fc Chakrabartil [20121 ). Another model of 
non-fluid bipolar outflows have been enthusiastically pur- 
sued by Becker and his collaborators which i nvolves Fermi 
acceleration of particles in isothermal sh ocks (|Le fc Becked 
120051 ; iBecker et a/.ll2008l ; IDas et aZ.ll2009l ). The advantage of 
shock in accretion model is that the presence or absence 
of shock, is good enough to broadly explain many varied 
aspects of black holes candidates such as spectral states, 
QPOs, jets an d outflows etc and the correlation between 
these aspects |Chakrabarti et all l2009h . Such corre l ations 
are also reported in observations (|Smith et aI1l200ll , |2002| ; 
iGallo et. al]|2003l ). 

There are other jet generation models too. Apart 
from ADIOS, magnetically driven outflows wer e also 
proposed by many authors ( Blandford fc Znaiek Il977l ; 
iBIandford fc Pavnefll982l ; IProgal 120051 ; lHawlev et al.ll2007l ). 



These magnetic bipolar outflows may be powered by the 
extraction of rotation energy of the black hole in the form 
of Poynting flux. However, recent observations finds weak 
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or no correlation of black hole spin with the jet formation 
around microq uasars and are conje ctured to be similar for 
AGNs as well (|Fender et gZ.I |2010| ). Bipolar outflows may 
even be powered by centrifugal or magnetic effect, and has 
been shown that for low angular momentum even weak mag- 
netic fields can produce equatorial inflow, bipolar outflow, 
polar funnel inflow and polar funnel outflow, and magnetic 
effect was identified as the main driver of such outflows. In- 
teresting as it may be, however, its connection with spectral 
and radiative state of black hole candidates is not well ex- 
plored. It is well known that the spectra of the black hole 
candidates extends to high energy domain (~ few MeVs), 
and one of the ways to generate such high energy non- 
thermal spectra is by shock acceleration of electrons, and 
has b een used to explain spectra from black hole candi- 
dates (|Chakrabarti fc Mandalll2006l ; iMandal fc Chakrabartil 
|200ST ) . Not only hydrodynamic calculations, even magneto- 
hydrodynami c investigations bears the possibility of shoc k 
in accretion l|Nishikawa et aZ.ll2005l ; iTakahashi et aZ.ll2006f ). 
Therefore we are looking for solutions of outflow generation 
which incorporates shocks, although presently investigating 
the effect of viscosity in formation of such outflows. 

Various simulations w it h Paczynski-Wiita potentia l 
(|Molteni et a/.lll994 Il996bl lal; lLanzafame et al.lll998l. l200ct) , 
as well as GRMHD simulations i Nishikawa et al\ 12003 ) 
showed the presence of post-shock bipolar outflows, which 
are far from isothermal approximation. Therefore, presently 
we relax the strict isothermality condition, and concen- 
trate on conservation of fluxes across the shock front. 
Theoretical frame work of thermally driven outflows has 
been done either for inviscid discs dChakrabartil Il999l; 
iDas et a?.ll200lj ; ISingh fc Chakrabartil I2012T I. or for viscous 
discs where the viscous stres s was assumed to be propor- 
tional to the total pressure (|Chattopadhvav fc Das! 120071 ; 
IDas fc Chattopadhvavl [20081 ). No theoretical investigations 
has thus far been made to study thermally driven outflows 
from discs where the viscous stress is proportional to the 
shear, a form of viscosit y which is probably more re alis- 
tic for black hole system ( Becker fc Subramaniumll2005l ). In 



Or 



this paper we solve viscous accretion disc equations for the 
contentious viscosity prescription, to compute the solution 
topologies of thermally driven outflows and the mass-outflow 
rates. We compute the energy angular-momentum parame- 
ter space for the accretion flows which will produce such out- 
flows, and its dependence on viscosity parameter. In numer- 
ical simulation with Shakura & Sunyaev viscosity prescrip- 
tion, shock location seems to increase with the increase o f 
viscosity parameter (|Lanzafame et aZ.ll998l ; lLee et qZ.l201lT ). 
while theoretical investigations with Chakrabarti-Molteni 
viscosity prescription showed the opposite phenomeno n 
IChattopadhvav fc DasfeoOTl ; IDas fc ChattopadhvavlkoOSD . 
We address this issue and attempt to remove any ambiguity. 
Since the post-shock disc could well be the elusive Compton 
cloud, while the Shakura-Sunyaev viscosity might be more 
physical since it satis fies proper inner b oundary conditions 
around black holes jBecker fc Lei 120031 ) . studies of shock 
driven outflows and its dependence on viscosity would throw 
light in understanding the radio X-ray correlation from X- 
ray binaries, especially, whether the jet becomes stronger or 
weaker with the increasing viscosity parameter, this would 
have an interesting connotation in interpreting observations. 
Since the post-shock disc can produce high energy photons, 
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Figure 1. Schematic diagram of accretion-ejection flow geometry. 
The black hole (BH) with its horizon is shown at O. OO 1 is the z 
axis, tcb = OM and rpw = OM' are spherical radial distances 
of the Centrifugal Barrier (CB) and Funnel Wall (FW). Moreover, 
= O'M', xcb = O' M are indicated. The two physical 
critical points x c i and x co of the accretion disc, and the shock 
location x s is also shown. The jet flows through CB and FW. 



then some part of the thermal energy gained through shocks 
will be dissipated as radiations. We study the issue of origin 
of such outflows in presence of dissipative shocks too. 

In the next section, we present the simplifying assump- 
tions and equations of motion. In section 3, we present the 
methodology of solution. In section 4, we present the so- 
lutions, and in the last section we present discussion and 
concluding remarks. 



2 MODEL ASSUMPTIONS AND EQUATIONS 
OF MOTION 

It has recently been stressed that black hole rotation plays 
no, or little part in generating or powering these jets 
(|Fender et al\\201(l ). So it is expected that plasma processes 
or fluid properties of the disc would generate jets. In this pa- 
per, we consider a non-rotating black hole and focus only on 
the fluid properties of accretion disc, which may be consid- 
ered to be responsible for jet generation. We have assumed 
the axis-symmetric disc-jet system to be in steady state. The 
space-time properties around a Schwarzschild black hole is 
described by the ps e udo-N ewtonian potential introduced by 
IPaczvriski fc Wiital <jl980h . The viscosity prescription in the 
disc is described by the Shakura-Sunyaev prescription. We 
ignore any cooling mechanism in order to focus on the effect 
of viscosity. The jets are tenuous and should have less dif- 
ferential rotation than the accretion disc, as a result the vis- 
cosity in jets can be ignored. To negate any resulting torque, 
the angular momentum at the jet base is assumed same as 
that of the local value of angular momentum of the disc. The 
accretion disc occupies the space on or about the equatorial 
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plane. But the jet flow geometry is described about the axis 
of symmetry. We first present the equations of motions of the 
accretion disc and the jet separately in the subsequent part 
of this section, and then in the next section, we describe the 
method to obtain the self consistent accretion-ejection solu- 
tion. In this paper we have used the geometric unit system 
where, 2G = M = c= 1 (M is the mass of the black hole, G 
is the Gravitational constant). Therefore, in this representa- 
tion the units of length, mass and time are the Schwarzschild 
radius or r g = 2GM/c 2 , M, and r g /c = 2GM/c i , respec- 
tively, consequently the unit of speed is c. 



Considering hydrostatic equilibrium in vertical direction, the 
local disc half height is obtained as: 



h(x) = t —ax 1 ^ 2 (x — 1) 
V 7 



The adiabatic sound speed is defined as 



(7) 



(8) 



where 7 is adiabatic index. The expression of the entropy- 
accretion rate is given by, 



2.1 Equations of motion for accretion 

The equations of motion for viscous accreting matter around 
the equatorial plane, in cylindrical coordinates [x, (j>, z) are 
given by, 

the radial momentum equation: 

1 \ 2 (x) 



du 1 dp 1 
dx p dx 2{x — l) 2 

The mass-accretion rate equation: 

M = 2-nT.ux, 







(1) 



(2) 



The mass accretion rate M is a constant of motion, except 
at the regions from where the mass may be lost into the jets. 
We present the exact form of the conservation of M later in 
the paper. The angular momentum distribution equation: 



dX(x) 
dx 



1 d(x 2 W x<t> ) 



Ex 



dx 







The entropy generation equation 

, da 



EuT- 



rf.r 



Q + -Q- 



(3) 



(4) 



The local variables u, a, p, p and A in the above equations 
are the radial bulk velocity, sound speed, isotropic pressure 
and specific angular momentum of the flow, respectively. 
Here, E = 2ph and W x <t> are the vertically integrated den - 
sity and the viscous stress tensor (|Matsumoto et aZlll984l ). 
Other quantities like the the entropy density, the local tem- 
perature, and the local half height of the disc are given by, 
s, T, and h, respectively. The local heat gained and lost by 
the flow are given by Q + — W 2 ^l(r\) and Q~ . 

The constant of motion of the flow is obtained, by inte- 
grating equation fTJ with the help of equations ([2] — 0, we 
find that the energy per unit mass is given by 



E 



2 2 
u a 

~2 + 7-1 



2x 2 



An A 



+ $0), 



(5) 



and is also called as the specific grand energy of the flow, and 
is conserved thro ughout the flow even in the presence of vis- 
cous dissipation (|Gu fc Lull2004 ). except across a dissipative 
shock (see, section 3.1.1). In Eq. flSJl, $(x) = -0.5/(x - 1) 
is the pseudo-Newtonian gravitational potential. 
The viscous stress is given by 

dn 

dx ' 



(6) 



where, r\ = pvh is the dynamic viscosity coefficient, v = 
aa 2 /(7$7fe) is the kinematic viscosity, a is the Shakura- 
Sunyaev viscosity parameter, and ilk are the local angular 
velocity and local Keplerian angular velocity, respectively. 



M{x) 



(9) 



If there is no viscosity i.e., a = 0, then M is constant except 
at the shock. At shock the entropy accretion rate will suf- 
fer discontinuous jump. The immediate pre-shock and post 
shock entropy-accretion rate denoted as M- and M+, are 
related as M+ > M- . But for < a < 1 or viscous flow, M 
varies continuously in the disc since viscosity dissipates and 
increases the entropy. If shock exists in viscous flow then, 
similar to the inviscid case M+ > M-. 

The gradient of the angular velocity can be obtained by 
integrating equation @ and also by utilizing equation 
and the expression of W X r/>, 



dM 
dx 



7«£4(A - An) 



(10) 



where, An is specific angular momentum at the horizon ob- 
tained by considering vanishing tor que at the event horizon 
l|Weinberdlf97ilBecker et aZ.ll2008h . Since A = x 2 ft, the ra- 
dial derivative of A is given by 



d\ _ 2 dQ 

— = 2xn + x — . 

dx dx 



(11) 



Moreover, fife denotes the Keplerian angular velocity and 
defined as 

1 



fi fc (a:) = 



2x(x- l) 2 ' 

The Keplerian specific angular momentum is defined as 



(12) 



Afe = fife x 



2{x - l) 2 



1/2 



(13) 



Manipulating equations |T] — with the help of equa- 
tions(0 — 113[) we obtain, 



du 
dx 



N 

15' 



where, 



N 



and 



2 (5a:-3)u (A 2 



7 + 1 2x(x - 1) 



7+1 



(14) 



w 2 Afe(A- A ) 2 



a 2 



7+1 

The gradient of the sound speed is, 



da (a ju \ du 
dx \u a J dx 



du (5x - 3)a 7(A-A 2 ) 



2x(x - l) 



(15) 



Therefore, the accretion disc problem in vertical equi- 
librium is solved by integrating Eqs. QUI 1141 US) . 
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2.1.1 Critical point conditions for accretion 

At large distances away from the horizon the inward velocity 
is very small and therefore the flow is subsonic, but, matter 
enters the black hole with the speed of light and therefore 
it is supersonic close to the horizon. Hence accreting mat- 
ter around black holes must be transonic, since it makes a 
transition from subsonic to supersonic. Therefore, at some 
location the denominator D of Eq. (|14[) , will go to zero, and 
hence the numerator iV goes to zero too. Such a location 
is called the sonic point or critical point. The critical point 
conditions are given as: 



M 2 = ^ 



7+1 



(5x c - 3)M C 2 



2x c {x c 
I 7 



a 3 c + 



(Ac — X 2 ) 



+7 



7+1 



M c A fe (A c - A ) 2 



(16) 



(17) 







where M c ,u c ,a c ,x c and A c are Mach number, the bulk 
velocity, the sound speed, the radial distance and the specific 
angular momentum at the critical point, respectively. 

The radial velocity gradient at the critical point is cal- 
culated by employing the l'Hospital rule. 



f dN/dx \ 
\dD/dx), 



(18) 



and by combining Eqs. (ll5l fc [T5)) we get, 



da\ _ f Oc _ 7^c\ ( du\ (5x e - 3)a c | 7(A 2 - A 2 . e ) 
dx ) c \u c a c J \dx J c 2x c (x c — 1) a c x| 

(19) 

So, the solution of Eqs. 1141 115|l . can only be ob- 
tained if we know the sonic point and its conditions (Eqs. 

IBS-EES. 



2.2 Equations of motion for outflows 

The flow geometry for accretion is described about the equa- 
torial plane, however, the jet or outflow geometry is about 
the axis of symmetry. If the outflow posses some angular 
momentum, then the outflow g eometry should be ho llow. 
Indeed, numerical simulations bv lMolteni et al] (|l996al ) sug- 
gests that the outflowing matter tends to emerge out be- 
tween two surfaces namely, the funnel wall(FW) and cen- 
trifugal barrier(CB). In Fig.l, the schematic diagram of the 
jet geometry is shown. The centrifugal barrier(CB) surface 
is defined as the pressure maxima surface and is expressed 
a.s 



x C b = [2X 2 rcB(rcB - l) 2 ] 4 



(20) 



where, rcB = \/ x 'cb + V~cb> spherical radius of CB. Here, 
xcb and ycB are the cylindrical radius and axial coordinate 
(i.e., height at res) of CB. We compute the jet geometry 
with respect to ycB i.e., %/fw — Vj = Vcb, where yFw and 
yj are the height of FW and the jet at rcB, respectively. 
The FW is obtained by considering null effective potential 
and is given by 



N2 (A 2 -2) + v /(A^-2)-4(l^y 2 J 

Xpw — A (AL) 



where, xfw is the cylindrical radius of FW. We define the 
cylindrical radius of the outflow 



xfw + xcb 



(22) 



The spherical radius of the jet is given by Tj = y^E 2 + y 2 



In 



3 1 »3 

Fig.l, OB(= r-j) defines the streamline (solid) of the outflow 
The total area function of the jet is defined as, 



.4 



2n(x 2 c 



x 2 FW ) 



\A + i dx il d Vi 



(23) 



where, the denominator is the projection effect of the jet 
streamline on its cross-section. The integrated radial mo- 
mentum equation for jet is given by, 



c 1 2 i 2 , ^3 
£j = + Ua 3 + 7^2 



2fo - 1) 



(24) 



where, £ j is the specific energy, Aj is the angular momentum 
of the jet, and n = 1/(7 — 1) is the polytropic index. The 
integrated continuity equation is, 



M out = pjVjA 



(25) 



and the entropy generation equation is integrated to obtain 
the polytropic equation (pj — Kjpj) of state for the jet. The 
entropy accretion rate for the jet is given by 



(26) 

In Eqs. (|24I25|I . the suffix 'j' indicates jet variables, where Vj, 
aj , and pj are the velocity, the sound speed and the density 
of the jet. 

Equations (|24I25I ) are differentiated with respect to 
r(~ vcb), to obtain, 



where, 



and, 



dr 



ilr , 



A 2 dxj 



2(rj 



dr 



dr 



aq dA 
~A~dr~' 



The critical point(rj c ) condition for jet is, 



2fo 



dr 



A 



dxj 
dr 



(27) 



(28) 



(29) 



Ac V dr 
(30) 



The outflow solution is obtained by integrating Eqs. 
(|27l) . with the help of Eq. (|30[) . The outflow equations can 
be determined uniquely if £j and Xj are known, because, the 
value of £j and Xj determines rj C . However, for consistent 
accretion-ejection solution, both accretion and jet part have 
to be solved simultaneously, a technique previously used by 
IChattopadhvav fc Dad (|2007f ). and has been refined in this 
paper. 



3 METHODOLOGY 

The accretion-ejection solutions are self consistently and si- 
multaneously solved, and we present the methodology to 
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find such solutions in this section. We start with the solution 
of the accretion disc. It has been mentioned before that Eqs. 
HUHUHS} can be integrated, if we know the sonic point x c . 
One of the long standing problem in accretion physics is to 
determine the sonic point of the flow in presence of a viscous 
stress of the form Eq. (0 , which keeps the angular momen- 
tum equation in a differential form [ Eg. (1 11 111 rather than 
an algebraic form (|Chakrabartilll996l ; ICu fc Lull2004l ). The 
problem is compounded by the fact that although quanti- 
ties on the horizon are known, but the coordinate singularity 
on the horizon makes it difficult to solve the equations by 
taking those as the starting values. The problem could be 
circumvented if the asymptotic values of the flow variables 

are known close to the horizon. 

W ith a clever use of conservation eauations lBecker fc Lei 
(|2003l ) found the asymptotic distribution of the specific an- 
gular momentum and the radial velocity close to the horizon, 
and they are, 



X(x) = A 



1 + 



2a f 2 



(2\ 1/2 {M 2 \~, N ^+<L 



(31) 



and 
u(x) = Uff(x) 



1 | 2Ex 2 -Ag-(7+l)/W 
x 2 u 2 ff (x) - (7 - l)f(x) 



1/2 



x — > 1 
(32) 



where the function f(x) is f(x) = -If^ [ 23 ,a^?~ Tg) ] 7+1 and 
the free fall velocity in the pseudo-Newtonian potential ge- 
ometry is given by 



(33) 
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Figure 2. The domain for multiple critical point MCP (dashed) 
and shock (solid) in E — A parameter space (e) is presented. The 
solution topologies or the plot of Mach number M with log(x) of 
the O type for E = 0.001, A = 1.5 (a); A type for E = 0.001, A = 
1.68 (b); W type fort E = 0.001, A = 1.75 (c); and I type for 
E = 0.005, A = 1.77(d), are presented. The solution type and 
their location in the parameter space is also indicated in (e). All 
the plots are for a = 0. The dotted lines in the panels named a, b, 
c, and d show all possible solutions while, the solid line with arrow 
heads show the actual accretion solution. The vertical long-short 
dashed line shows the location of the sonic points. 



3.1 Integration procedure of accretion solution 

A position very close to the horizon X — Xin — 1.01 is chosen. 
For given values of the parameters a, E, Xq at the horizon 
and Mm at x in , Eqs. (f3Tj) , ((32} and (f33)l can be used to de- 
termine the asymptotic values of the fluid variables close to 
the horizon, and which can be used as the initial values for 
the integration. With these values, Eqs. (|11II14|[T5)) are inte- 
grated outwards, simultaneously checking for the location of 
the sonic point by using Eqs. fj 161 — I19fl . The sonic point x c is 
determined iteratively, and once it is found, the sonic point 
conditions are used to integrate the equations of motions, 
outwards. It is well known that matter with angular mo- 
mentum may posses multiple critical points (MCP). But the 
flow may pass through various sonic points only if a shock 
is present, in other words, the existence of MCP is manda- 
tory for the formation of shock in black-hole accretion. Only 
a supersonic flow can undergo shock transition and so exis- 
tence of one sonic point at larger distance (i. e.,x co ) from the 
horizon is warranted. The post shock flow is subsonic. How- 
ever, the inner boundary condition of black hole accretion 
is supersonic, and hence the subsonic post-shock flow has to 
pass through another sonic point (x c i), before it dives into 
the black hole. In other words, shock in black hole accretion 
may exist only if MCP exists. It is to be remembered though, 
that there is no smooth transition between various branches 



of the solution passing through different sonic points. This is 
true for both inviscid and adiabatic flows as well as viscous 
and non-adiabatic flow. However, if the flow is inviscid, for a 
given value of E and Ao all possible sonic points are known 
apriori. Incase of viscous flow which is following a viscosity 
prescription of the form Eq. ([6}, the existence of multiple 
sonic points can only be ascertained only if there is a shock. 



3.1.1 Shock equations 

The Rankine Hugoniot (RH) shock conditions are obtained 
by conservation of the mass, momentum and energy flux 
across the shock. In presence of mass loss and energy loss, 
the shock conditions are given by, the modified mass con- 
servation, 

M+ = M- - M out = AL{1-R 7h ), (34) 

This equation effectively divides the accretion rate of pre- 
shock accretion disc (M_) into two channels, namely the 
post-shock accretion disc (represented by M+) and the jet 
(represented by M out ). The modified momentum conserva- 
tion, 

W+ + Y.+U+ = W- + (35) 



© 2012 RAS, MNRAS OOO.mflBI 



Mass outflow from viscous discs 7 



and the third shock condition is the modified energy conser- 
vation 

E+ = E-- AE, (36) 
where, Rm is the relative mass outflow rate given by, 

M out 



R-ri 



M- 



(37) 



Here, subscripts minus(-) and plus(+) denote the quan- 
tities before and after the shock. W is the vertically 
integrated pressure. In absence of massloss (R,h = 0) 
or dissipation (i.e.,AE = 0), Eqs. (j34] |35] [36} reduce 
to the standard RH shock conditions. Since the dissi- 
pation is expected at the shock location, it is assumed 
that energy dissipation t akes place mostly through the 
thermal Comptonization (|Chakrabarti fe Titarchukl 1 19951 ; 
iMandal fe Chakrabartill2010l ) and is likely to be very im- 
portant within a distance dx inside the shock where the op- 
tical depth is around unity. So this energy dissipation in the 
post shock flow reduces the temperature of the flow and the 
loss of energy is proportional to the temperature difference 
between the post-shock and the pre-shock flows, i.e. 



AE = f e n(a 2 + - a_) 



(38) 



where f e is the fraction of the difference in thermal energy 
dissipation across the shock transition and n is the poly- 
tropic index. We use f e as a parameter, in presence of de- 
tailed radiative processes f e can be self-consistently deter- 
mined. Since shock width is infinitesimally thin, so we as- 
sume that dQ/dx is continuous across the shock. The angu- 
lar momentum jump condition is calculated by considering 
the conservation of angular momentum flux, and is given by 



A_ = A+ + C a 



where, C s h = —u+(X^ 
written as, 



Ao)/<4- 



Equation 



(39) 



can be re- 



(A- - A ) 



<2-M+(A+ — Aq) 



(40) 



Since at shock, a+ > et_ and u_ > u+, therefore, A+ > A_. 

Using shock condition equations (|34[) and (j35j), the pre- 
shock sound speed and bulk velocity can be written as, 



2 , 2 
a_ = Ki«- — "/u_ 



(41) 

where, k\ = {a\ + ^u 2 + )/{fu+), f = 1/(1 — Rm)- Now sub- 
stituting for a_ and A_ in equation (|36|l . we find a quadratic 
equation of u~ as, 



where, 



C*2 



C 2 u_ + Ciu_ + Co = 



1 . 7 2 C? h 



(42) 



"7(1 + fe) 



2x1 



Ci = 



nil + f e )ki 1 — 



Co 



A c _ <y_ 
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Figure 3. Variation of M with log(x) for the accretion solutions 
with different viscosity parameter a. In a, c. i, m we present in- 
viscid solutions corresponding to the O, A, W, I type of solutions 
from Figs. 2a-d. Towards right E, \q is kept the same but a is in- 
creased. The flow parameters for which these plots are generated 
E = 0.001, A = 1.5 (a, b, c, d); E = 0.001, A = 1.68 (e, f, g, h); 
E = 0.001, A = 1.75 (i, j, k, 1) and E = 0.005, A = 1.75 (m, n, 
o, p). The viscosity parameter a is mentioned on the figure. The 
vertical long-short dashed line shows the location of sonic points. 



C = Aj 



' Csh 



U+ 



-hi 



and 



k 2 = E+ - $(x- 3 ) + f e na + . 

In terms of the shock quantities, the mass outflow rate 
is given by 

Rvj(x s )A(x s ) 



Rrh = M out /M_ = 



4tt 



3/2/ tn 

x a (x s — l)a+u- 



= RVQ, (43) 



where, R — E+ /£_ is the compression ratio across the 
shock, V = Vj(x 3 )/u- is the ratio of jet base velocity 
and the pre-shock velocity of the disc (it_), and Q — 
A{x a )/[A-R{2/ 1 ) 1 ' 2 x 3 J 2 {x s -l)a + ]is the ratio of the jet cross- 
sectional area at r — x s and the post-shock accretion disc 
cross-sectional area. 

3.2 Accretion-Ejection Solution 

The accretion-ejection is computed self-consistently. We 
have set 7 = 1.4 and Xi„ = 1.01 throughout the paper. The 
method to find the accretion-ejection solution is as follows, 

(i) First we assume that R,h = 0. With chosen values 
of Ao, E, a and f e , we follow the procedure described in 
section 3.1, i.e., determine the inner sonic point iteratively 
and integrate outwards. Equations (|39l 1411 142[) are checked 
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Figure 4. Comparison of shock free (dotted) and shocked solu- 
tion (solid). Both the solution starts with the same outer bound- 
ary condition Xi n j = 5117, Xi n j = Ax(xinj) = 50.6, E = 10 — 3 . 
The viscosity parameter for the shock free solution is a = 0.02 
(dotted) and that for the shocked one a = 0.15 (solid). Various 
flow variables are plotted are M (a), a (b), u (c), A (d), E (e) 
and Q (f) as a function of x. The sonic points are marked with 
open circles (a). The shock free solution has only one x c i, the 
shocked solution has both x c i and x co and the shock location 
is at x s = 49.68. Dashed plots of A^- and Qfc in (d) and (f) 
are Keplerian angular momentum and angular velocity drawn for 
comparison. 



to calculate the pre-shock quantities. We find out the outer 
sonic point (x co ) iteratively, from the pre-shock quantities. 
Two possibilities may arise, either the flow passes through 
only one sonic point and gives a smooth solution, or, finds a 
stable shock solution and passed through two sonic points. 
The location of the jump for which x co exists is the virtual 
shock location (x s ). 

(ii) Once x B is found out, we assign £j = £+ and Xj = A+ 
and solve jet Eqs. (|27l fc[30 )l and compute the corresponding 

Rrh- 

(iii) We put this value of Rm into Eqs. (|39I42[) to find a 
new shock location . 

(iv) Steps (ii) & (iii) are repeated till the temporary shock 
location x s converges to the actual shock location (x„). The 
converged shock solution therefore gives the actual jet solu- 
tion too. We find that the Mj > M- and M+ > M-. Since 
matter prefers higher entropy solutions, therefore, the post- 
shock fluid would prefer both the channels, one through the 
x C i onto the black hole, and the other through r JC and out 
of the disc-jet system. In other words, a shock in accretion 
would generate a bipolar outflow from the post-shock region 
too. 

The outer boundary is chosen as Xinj = 10 4 7" 9 , or the dis- 
tance at which X(xinj) = \ic{xinj), which ever is shorter. 



Figure 5. Variation of A with log(x) of shocked accretion flow 
(a, b, c), for parameters E = 0.0001, & A;, = 1.68. Each curve 
represents a = 0.0 (solid), a = 0.0075(dotted) and a = 0.015 
(dashed). First three panels are for x^ = 2220r 9 (a), x^ = 3660r 9 
(b), and Xf, = 5230r 9 (c). Vertical solid line and dash-dotted line 
show the shock location for a = 0.0 curve and outer boundary 
location (xj), respectively. The variation of limiting x;, with E 
for A b = 1.68 (d), and A;, = 1.65 (e). Domain 1 represents all x b 
for which x s decrease with the increase with a, but for any Xj, in 
2, x a will increase with the increase of a. 

4 SOLUTIONS 

We present every possible way matter can dive into the black 
hole. In this section, we start with accretion solutions with- 
out considering massloss or dissipation at the shock front 
and study the effect of viscosity on accretion solution. Then 
we present the accretion-ejection solution and study how the 
viscosity can affect the mass outflow rates. And finally we 
present accretion-ejection solutions in presence of dissipa- 
tive shocks and show the effect of both the viscosity and 
dissipation at the shock front. 

4.1 All possible accretion solutions in the 
advective regime 

Since black hole accretion is necessarily transonic, at first, we 
present the simplest rotating transonic solutions ie.,inviscid 
global solutions (global solutions which connect horizon and 
large distances) in Figs.[2^-d. The inviscid solutions are pre- 
sented as an attempt to recall the simplest accretion case 
(|Chakrabarti| [l989; iDas et all 120011 ). For inviscid flow the 
constants of motions are E, A (= Ao), and if outflows are 
not present then, M is also a constant of motion. Moreover, 
in absence of viscosity, it straightaway follows from Eq. ([SJ 
that E = £ = u 2 /2 + na 2 + \ 2 /(2x 2 ) - 0.5/(a: - 1), where £ 
is the Bernoulli parameter or the specific energy of the flow. 
In Figs. (2K _ d, accretion solutions in terms of the Mach num- 
ber M (= u/a) distribution are presented, and in Fig. [2]3, 
the E — Aq parameter space for multiple sonic point and the 
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Figure 6. E — Ao parameter space for shock, for various viscosity 
parameters a = 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3 marked on the 
figure. No energy dissipation (/ e = 0.0) at the shock and no mass 
loss (/ = 1.0) is assumed. 

shock is presented too. Depending on E & Ao of the flow, the 
solutions are also different. If the Ao is low, there is only one 
outer sonic point x co , and solution type is O-type or Bondi 
type (E — 0.001, Ao = 1.5, Fig. As Ao is increased, the 
number of physical sonic points increases to two and the ac- 
cretion flow which becomes supersonic through x co can enter 
the black hole through x C i if a shock condition is satisfied. 
Although the shock free solution is possible but in this part 
of the parameter space a shocked solution will be preferred 
because a shocked solution is of higher entropy (or in other 
words of higher M). Such a class of solution is called A- type 
(E = 0.001, Ao = 1.68, Fig. HJa). For even higher A only 
one sonic point is possible (E — 0.001, Ao = 1.75 for Type 
W; and E = 0.005, A = 1.77 for I shown in Figs Hi and 
d), and the solutions are wholly subsonic till x c i and then 
dives on to blackhole supersonically. W type solutions are 
different than I type, in the sense, W type is still within 
the MCP domain while I type is not. Moreover, I type is 
a smooth monotonic solution, although W is smooth and 
shock free but is not monotonic and has an extremum at 
around x co . The parameter space E — Ao, bounded by solid 
line (mno) shows the RH shock parameter space, while the 
dotted one (PQR) shows the MCP domain (Fig. 2e). It is to 
be noted, that in the inviscid limit, accretion is only possible 
if A < A A '. 

Figures [3^, [3^, [3j, [3jxi, represent the invis- 
cid solutions corresponding to the O type solutions 
(E, Ao = 0.001, 1.5), A type (E, A = 0.001, 1.68), W- 
type(S, Ao = 0.001, 1.8), and I-type (E, A = 0.006, 1.8), 
and are also depicted in Figs. 2a-d. Keeping E & Ao same, 



we increase the viscosity parameter in the right direction. 
Figures [3^,-d, has same E & Ao, but progressively increasing 
a = 0.06 (Fig. EJd), 0.068 (Fig. W) and 0.1 (Fig. [3H). 
Similarly for Figs. [3^-h, E — Ao is same but different a, 
so is the case for Figs. [3J-1, and Figs. (3jn-p. Interestingly, 
the viscous I type is in principle the much vaunted ADAF 
type solution presented in Figs 3d, 3h, 3k-l, and 3n-p. The 
shock-free solution is characterized by monotonic spatial 
distribution of flow variables, and wholly subsonic except 
very close to the horizon which are essentially th e viscous 
I typ e solutions and has also been shown by iLu et al\ 
\ 19991 ) ; lBecker et al\ (|2008l ); iDas et al\ (|2009T l. It is evident 
from Figs. [3]i-p, that the effect of viscosity is to create 
additional sonic points in some part of the parameter 
space, opening up of closed topologies, and might trigger 
shock formation where there was no shock, while removing 
both shock and multiple critical points in other regime of 
the parameter space. All of this is achieved by removing 
angular momentum outwards while increasing the entropy 
inwards. In this connection one may find two kind of 
critical viscosity parameters in the advective domain. If the 
inviscid solution is O type (Fig. 3a), then there would be a 
lower bound of critical viscosity a c i which would transport 
angular momentum in a manner that would trigger the 
standing shock. And there would be another upper bound 
of viscosity parameter q cu that would quench the standing 
shock. While, if the inviscid solution has a shock to start 
with (Fig. 3e), then there could only be a cu . For the case 
presented in Figs. 3a-d, a c \ — 0.0465 and q cu = 0.065. And 
for the case presented in Figs. 3e-h, a cu — 0.0126. 

So far, we have compared solutions with different vis- 
cosity but same internal boundary condition (E, Ao). Let 
us compare two solutions with same outer boundary con- 
dition. In Figs. 0^-e, we compare two solutions starting at 
the same outer boundary Xi n j = 5117, same grand energy 
E = 10~ 3 and \i n j — \K{xinj) = 50.6. Both the solution 
starts with the same entropy ~ Mmj = 2.72 x 10~ J . For 
a = 0.15 (solid) the accretion solution has a standing shock 
at x s — 49.68 (the vertical jump in solid), while the shock 
free solution is for a = 0.02 (dotted). The flow variables 
plotted are M (Fig. H*), a (Fig. [4p), u (Fig. gt), A (Fig. 
HJi), and £ (Fig. 0^). Incase of inviscid flow E = £ is con- 
stant, for viscous flow E is still a constant of motion but 
£, the specific energy, varies with x as is shown in Fig. [4^. 
The inner part of the shocked solution is faster, hotter, and 
of lesser angular momentum. The inner boundary for the 
shocked flow is E = 10~ 3 , Ao = 1.29, while that for the 
shock free solution is E = 10~ 3 , Ao = 1.7. The jump in A at 
the shock follows Eg. 1401 Interestingly, the post shock flow is 
of higher entropy [A4(a; ~ 1) = 7.49 x 10 -5 ] than the inner 
part of the shock free solution [M(x ~ 1) = 7.12 x 10~ 5 ]. 
Although it may seem contradictory that a shock free solu- 
tion exists for lower a, while shocked solution appears for 
higher a, however, it is to be remembered that with such 
high Xinj, there would be no accretion solution for a — 0. 
One has to have a certain non-zero viscosity to even have 
a global (that connects horizon and outer edge) solution. 
Infact, for this particular case we have identified a limiting 
viscosity parameter Qi = 0.00747, such that advective global 
solutions are possible for any viscosity parameter a ^ ai. 
And for flows starting with such high initial A, if the vis- 
cosity parameter a is small, then the A distribution will be 
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Figure 7. Mach number distribution for self-consistent accretion- 
ejection solution (a & b). The disc solution (a), and the jet (b) 
are plotted for accretion flow parameters E = 0.001, Ao = 1.542 
and a = 0.05. The shock in accretion is shown by the vertical 
jump at x s = 56.236 (a). The density distribution pj (c) and 
the temperature Tj distribution of the jets are also plotted (d). 
The sonic points arc marked by open circles (a & b). The dotted 
vertical line shows the position of the shock if R^ = 0. The jet 



flow parameters £j = £4 



1.084 X 10~ 3 , A, = A 4 



1.699. 



higher. Therefore, for such high A radial velocities will not 
be high enough to become supersonic and form a shock. 
So one needs higher a a c i = 0.1496 to reduce the an- 
gular momentum to the extent that may produce standing 
accretion shock. And if a is increased beyond another crit- 
ical value a cu = 0.1555, steady shock is not found. Hence 
one can identify three critical a's for accretion flows starting 
with the same outer boundary condition, where the angular 
momentum at the outer boundary is local Keplerian value. 
Interestingly, a prescription originally in voked to generate 
Keplerian disc |Shakura fc Sunvaevlll973h . can produce sub- 
Keplerian accretion flow with or without shock even if angu- 
lar momentum is Keplerian at the outer boundary. This is 
not surprising since, the gravitational energy released by the 
infalling matter, would be converted to kinetic energy and 
thermal energy (by compression and viscous dissipation). 
If the gravitational energy is converted to thermal energy 
and only the rotation part of kinetic energy, and also if the 
thermal energy gained by viscous dissipation is efficiently 
radiated away, then one would produce Keplerian disc so- 
lutions. The \k (Fig- StO and fix (Fig. U) are presented 
for comparison. In the present paper, the radiative processes 
have been ignored which produces hot flow, but the advec- 
tion terms have not been ignored. Therefore, sub-Keplerian 
flow with significant advection are obtained. 

The effect of a on shock location x 3 is an interesting is- 
sue. Numerical simulations show that, for fixed outer bound- 
ary conditio n, x a expands to larger distances with the in- 



crease of a (ILanzafame et al.\ 119981 ; Lee et al.\\2011u . while 
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analytically IChattopadhvav fc Daa (|2007T ) showed that for 



Figure 8. Various shock and jet quantities like R^ (a), x s (b), 
R (c), V (d), Q (e), and M (f) arc plotted with E for a = 0.05 
and Aq = 1.54. 



the same outer boundary condition, x„ shifts closer to the 
horizon wit h the increase in a. Although, the viscosity pre- 
scription of I Chattopadhyay &: Dasl (j2007n and the simula- 
tions are different, namely the former chose the stress to be 
proportional to total pressure, while in simulations the stress 
is proportional to the shear, still viscosity reduces angular 
momentum, and we know for lower angular momentum if the 
shock forms, it should form closer to the black hole! Since 
the viscosity prescription in this paper is similar to the sim- 
ulations, we should be able to answer the dichotomy. So a 
concrete question may arise, if viscosity is increased, does 
the Xs expands to larger distances or, contracts to a posi- 
tion closer to the horizon? If the viscosity acts in a way such 
that the A+ (immediate post-shock A) is less than its inviscid 
value at the shock then x s will move closer to the horizon. 
However, such simple minded reasoning may fail, if shocks 
exists, then the post shock flow being hotter would trans- 
port more efficiently than the immediate pre-shock flow. So 
although, A_ < \(xi n j), it is not necessary A+ will be less 
than X(xi„j). We scourged the parameter space to find the 
answer, and in the following we present the explanation. 

Let Xb be the distance at which the A distribution of the 
viscous solution is coincident with the A value of the inviscid 
solution, let us further assign A;, = X(xb). It is to be remem- 
bered that Ao = Xb for the inviscid (a = 0) solution, but 
Ao < At for viscous solution because viscosity reduces the 
angular momentum. In all the simulations done on viscous 
flow in the advective regime referred in this paper, one starts 
with an inviscid solution (since analytical solutions are read- 
ily available), and then the viscosity is turned on, keeping 
the values at the outer boundary fixed. In other words, it is 
this Xb that is called the outer boundary in numerical simu- 
lations, and generally, Xb ^ few x 100 since it is computation- 
ally expensive to simulate a large domain from just outside 
the horizon to few x 1000 r g and still retain required resolu- 
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tion to achieve intricate structures in the accretion disc. In 
Figs. [5^,-c, we compare the X(x) of shocked accretion flows 
starting with the same E = 10 -4 , & Xb = 1.68 for vari- 
ous viscosity parameters such as, a = (solid), a = 0.0075 
(dotted), a = 0.015, but for different points of coincidence, 
e.g.,x b = 2220 (5a), x b = 3660 (5b), and x b = 5230 (5c). The 
vertical solid line is the location of the shock for the inviscid 
flow. Although, we match the At, at xt of the viscous and in- 



viscid solutions, we still integrate outwards upto Xi 



10 4 



Since Xi„j >> 1, this distance may be considered as the size 
of the disc. In Figs. [5ji-c, we show that depending upon the 
choice of Xb, the outer boundary conditions can be remark- 
ably different. If Xb is short (Fig. [5^0, then flow with higher 
a will be able to match Xb at the same Xb, only if one starts 
with much higher \i n j, in other words the gradients will be 
stiffer. Consequently, the increase of a will create higher A+, 
which will result in the increase of x B . In case Xb is large (Fig. 
[S}i), the gradients are smoother and the resulting A;„j will 
be approximately similar for any value of a. Hence as a is 
increased, A+ would decrease and consequently x s will de- 
crease too. In Fig. [5^, the parameters are A; n j = A b = 1.68 
for inviscid or a = (solid), Xi„j = 1.769 for a = 0.0075 
(dotted) and Xi n j = 1.858 for a = 0.015 (dashed). Since 
Xb = 2220 is short, the gradients are steeper, and as dis- 
cussed above, in such cases x s increases with increasing a. 
While in Fig. X inj = 1.744 for a = 0.0075 (dotted) and 
Xinj = 1.807 for a — 0.015 (dashed), where x b = 5230 is 
larger, the shock x s decreases with increasing a. Since for 
shorter values of Xb, x s increases with the increase of a, and 
for longer Xb, x s decreases with the increase of a, so a lim- 
iting Xb should exist for which x s will neither increase or 
decrease with the increase of a. In Fig. [S]], we show that for 
Xb = 3660, the shock neither increase or decrease with the 
increase of a. 

In Figs. [5ji & e, we plot the limiting Xb as a function of 
E for Ao = 1.68 (d) and Ao = 1.65 (e). The domain name 
'1' corresponds to any Xb at which the x s will decrease and 
'2' signifies the domain where at any Xb, x a increases with 
a. Hence Fig. [S^, lies in domain 2, Fig. on the curve, 
and Fig.[5ji on domain 1 of Fig.[5ji. Incidentally, if the outer 
boundary condition for all the advective solutions start with 
the Keplerian angular momentum, then the shock location 
x s decreases with the increase of a. Since the numerical sim- 
ulations are usually performed with a shorter radial extent 
(xb ~ few x 10 — few x 100), this is similar to the case Xb 
lying in the domain 2. As a result, earlier simulations of ad- 
vective accretion flows have reported the increase of x s with 
the increase of a. Hence we may conclude that the incre- 
ment of x 3 with the increase of a is an artifact of faulty 
assignment of outer boundary condition in the simulations. 

In Fig. H3 the E — Ao parameter space for standing ac- 
cretion shock has been plotted for various viscosity param- 
eters, a = 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3. Since viscos- 
ity will in general reduce the angular momentum along the 
flow, Ao should decrease with the increase of a. As a result, 
the shock parameter space shift to the lower end of the Ao 
scale. One may compare the RH sho ck parameter sp ace with 
that of the isothermal shock space (|Das et ffll.||200g| ). For all 
possible boundary conditions, RH shocks may be obtained 
upto a = 0.3. For values outside the bounded regions there 
are no standing shocks, although transient and oscillating 
shocks may exist. It is to be remembered though, the pa- 
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Figure 9. E — Ao parameter space of the standing adiabatic 
shock with Shakura-Sunyacv viscosity parameter a = 0.0, 0.05, 
0.1, 0.15, 0.2 marked on the figure, and in presence of massloss. 



rameter space shown here corresponds to inner boundary. 
In the outer boundary A might be much higher for viscous 
flow as has been shown in Figs. [4^i-e. 



4.2 Inflow-Outflow solutions 

In the previous sub-section we discussed all possible solu- 
tions that one can have in presence of viscosity prescription 
given by Eq. ([6]), and non-dissipative or RH shocks. In this 
section we present the self consistent inflow-outflow solu- 
tions. In presence of massloss, the mass conservation equa- 
tion across the shock will be modified in the form of Eq. 
I|34|) . All the steps mentioned in section 3.2 are followed to 
compute the self-consistent inflow-outflow solution. In Fig. 
[7^-d, we present one case of global accretion-ejection solu- 
tion. In Fig. [T^l, the accretion solution in terms of M(x) 
or the Mach number is presented, in Fig. 03, the jet Mach 
number Mj is plotted with the height j/j of the jet from the 
equatorial plane of the disc, where the radial coordinate of 

the jet is rj = 



x? + Uj . The sonic points are marked with 
open circles. The inner boundary for the accretion solution is 
E — 0.001, Ao = 1.542. The specific energy and the angular 
momentum at the shock are £+ = 1.084 x 10" 3 , A+ = 1.699. 
The collimation parameter of the jet (see, Chattopadhyay 
2005) at its sonic point is Xj c /yj c ~ 0.24, and hence the 
spread is quite moderate. The dotted vertical line is the po- 
sition of the shock when i?„ = 0, and the solid vertical line 
is the position of the shock after the mass-outflow rate has 
been computed. Clearly, since excess thermal gradient force 
in the post-shock disc drives bipolar outflows, it reduces the 
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Figure 10. Comparison of accretion solutions with non- 
dissipativc shock i.e.,f e = (solid), and dissipativc shock f e = 
0.1 (dotted). The inner boundary condition for f e = arc 
E, A = 7.373 X 10~ 4 , 1.527 (solid), and for /,. = 0.1 arc 
E + , A = 2.5 X 10~ 4 , 1.527 (dotted). In case of the dotted curve 
E- = 7.373 x 10~ 4 . For / e = 0, x 3 = 20.15 , & R, h = 0.1133 
(solid), and for f e = 0.1, x s = 18.04 & R^ = 0.0755 (dashed). 
Various flow variables are M (a), u (b), A (c), £ (d), a (e), and 

M (f). 



pressure, and hence to maintain the total pressure balance 
across the shock, the shock front moves closer to the horizon. 
The relative mass outflow rate computed for the particular 
case depicted in Figs. [T^i-d, is Rm = 0.1044. In Fig. 0:, we 
plot the density pj of the jet derived for M = O.lMEdd for a 
black hole of M — IOMq . In Fig. [TJi, we plot the tempera- 
ture Tj of the jet. As is expected the jet is hot near its base 
but falls to low temperatures at larger distances away, while 
the density also falls as is expected of a transonic jet. 

Figures [8k,-f, are plotted for the fixed values of a = 0.05 
and Ao = 1.54. Various quantities Rm (a), x s (b), 7? (c), 
V (d), Q (e), and M (f) are plotted with E, which means, 
we are studying how Rm change with the change of inner 
boundary condition. Generally the Rm is few percent of the 
injected accretion rate but it can go upto more than 10% of 
M- for high E. Although in this plot, our main interest is to 
quantify the relative mass outflow rate Rm with E for given 
value of a and Ao, we have plotted all the essential quantities 
that would contribute in driving a part of the post-shock 
matter as outflows. The shock location x a increases with 
the increase of E, and for increasing x s , R should decrease. 
Since Rrh is combination of R, V and Q, even for falling 
R and Q, the mass outflow rate increases since it is being 
compensated by the increase of V. Interestingly, none of the 
quantities show a monotonic variation. It is important to 
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Figure 11. Relative mass outflow rate R^ as a function of (a) 
E— for fixed values of a = 0.02, Ao = 1.65; (b) Ao for fixed values 
of = 0.001 and a = 0.02; (c) a for fixed values of = 0.001 



and Xi„ 



70.6 at Xi. 



10 4 ; where all the plots are for f e 







(solid), f e = 0.05 (dotted) and f e = 0.1 (dashed), (d) R^ as a 
function of f e for fixed values of = 0.001 and Ao = 1.65 and 
each curves are for a = (solid), a = 0.01 (dotted) and a = 0.02 
(dashed). 



note, all the three parameters R, V and Q represent the 
jet to disc connection and not the actual driving. The real 
drivers are however, the post-shock specific energy £ + , and 
the jet base cross-section ^4 S (= A(x s )). Higher £+ means 
hotter flow at the jet base, and therefore the thermal driving 
will be more. This is complemented by the cross-section As 
of the jet base. Higher £+ would drive more matter into the 
jet channel but will be limited by the cross-sectional area. 
In Fig. [Hf, we plot the pre-shock entropy-accretion rate M- 
(solid), the post-shock entropy-accretion rate M+ (dashed) 
and the jet entropy-accretion rate Mj (dotted), as a function 
of E, and it is quite evident that the pre-shock entropy is less 
than both Mj and M+. Moreover, when the value of Mj 
is high, Rm is found to be high too, which vindicates that 
matter would prefer to flow through channels with higher 
entropy. 

From Eq. (JM}, we know that M + < M_ if R^, 0. 
Therefore, the post-shock pressure would be reduced. This 
would cause x s to dec r ease as shown in Fig. [7^, [also see 
IChattopadhvav fe Dasl (|2007l )]. Furthermore, the massloss 
from the post-shock flow and consequent reduction in pres- 
sure would also modify shock parameter space. In Fig. [9] the 
bounded regions in the E — Ao space, represents the shock 
parameter space for a — 0, 0.05, 0.1, 0.15 and 0.2, marked 
on the figure. Compared to the shock parameter space in ab- 
sence of massloss (Fig. , the parameter space in presence 
of massloss gets reduced and beyond a — 0.2 the stand- 
ing shock seems to vanish. Moreover, the shock parameter 
space shows that the standing shock do not seem to exist 
for very low E. Non-existence of standing shock ofcourse do 
not imply non-existence of non-steady or oscillatory shocks. 
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lently, E-, Ao, a, f e 



E-, Xi n j, a, f e . If f e = then 



E — E+ = E- , therefore any of the following two sets would 



suffice E, Ao, a, or, E, \i r , 



Since the outflows are 
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Figure 12. Relative mass outflow rate R m (solid), x s (dashed), 
£-l_ (dotted), and As (long-dashed) as a function of a for f e = 
(a) and f e = 0.1 (b) for fixed values of E- = 0.001 and \i n j = 
10 4 . 



4.3 Massloss from the dissipative shocks 



Chakrabarti fe Titarchukl (|l995l ) and later 

Mandal fe Chakrabartil (|2010h had shown that the post 



shock hot flows can produce the high energy photons easily 
by inverse- Comptonizing the soft photons and reproduced 
the observed spectra from variety of objects. If indeed 
the observed spectra and luminosity can be reproduced 
from the post-shock flow, then the grand energy will not 
be conserved across the shock and will be given by Eq. 
(|36[) , and the dissipated energy can be radiated away. In 
Fig. HOb-f. various flow variables M (a), u (b), A (c), £ 
(d), a (e), and M (f) are plotted for f e = (solid) and 
f e = 0.1 (dotted). For f e = (solid), the solutions are 
plotted for flow parameters E, A = 7.373 x 1CT 4 , 1.527. 
For f e = 0, we have E = E+ — E-. However, for f e > 0, 
E+ < E-. So for f e = 0.1 (dotted), the inner boundary 
is represented by E+, Ao = 2.5 x 10 -4 , 1.527, while the 
pre-shock E- — 7.373 x 10~ 4 . The post-shock specific 
energy or Bernoulli parameter, temperature, angular 
momentum, and entropy of the solution with dissipative 
shock (dotted) are lower than those corresponding to the 
non-dissipative shock (solid). Since a part of the thermal 
energy gained through shock is spent in powering jets and 
to produce radiation, the shock front moves closer to the 
black hole. The relative mass outflow rate or R m is lower 
for dissipative shocks, because part of the thermal energy 
of the post-shock disc is lost through radiation. However, 
the thermal energy lost as radiation can still contribute to 
j et power if those photons deposit momentum onto the jets 
(IChattopadhvav fe Chakrabartill2002l ; lchattopadhvav et all 
12004 IChattopadhvavll2005l ). 

If fe =/= and shocks are present, then the accretion 
solutions are obtained either for E+, Ao, a, f e , or equiva- 



launched from the post-shock flow, the mass outflow rates 
are computed only from the shocked accretion flows, and 
we are interested to find the dependence of R m on each of 
the above mentioned parameters. In Fig. Illfa . we plot R m 
as a function of E- when the other parameters are fixed at 
a = 0.02, Ao = 1.65 where each of the curve are for f e = 
(solid), fe = 0.05 (dotted) and f e = 0.1 (dashed). With 
the increase of E-, the post-shock specific energy £+ in- 
creases which drives more matter as outflow, and hence R m 
increases. However, for any given E-, R m decrease with the 
increase of f e , i.e., the mass outflow decreases with the in- 
crease of the thermal energy dissipation. Although for very 
high E-, the separation of the curves decreases, this shows 
that flows starting with higher energy would have enough 
thermal energy to drive significant outflows even for f e ^ 0. 
Now, let us fix E- but vary Ao in Fig. II 1 b . R m is plot- 
ted with Ao for f e = (solid), f B = 0.05 (dotted) and 
f e — 0.1 (dashed). The fixed parameters are E- = 0.001, 
and a = 0.02, and clearly Rm is not a monotonic function 
of Ao. Increasing Ao would increase A+ and therefore increase 
x 3 but would decrease £+. The increase of x s , increases the 
jet base cross-section. The competition between £+ and A 
causes a dip in Rm for moderate values of Ao- 

However, R m decreases with the increase of f e . In Fig. 
1 1 lb . we fix the outer boundary condition, namely, E- — 

0. 001, \i„j = 70.6 at Xi n j = 10 4 , and R m is plotted as a 
function of a, where each of the curves are for f e = (solid), 
fe — 0.05 (dotted) and f e — 0.1 (dashed). Since for flows 
starting with such high Xi n j there would be no accretion 
shock solution for low a, therefore one observes outflows only 
beyond a critical value of a. We know (see Fig. [5fi) that for 
flow starting with the same outer boundary condition, the 
shock location decreases with the increasing a. Decrease in 
x s would mean decrease in As- Since x s decreases, it means 
the post-shock energy increases due to viscous dissipation 

1. e.,£+ increases. Increase in £+ would drive more matter 
into the outflow channel. Therefore, the decrease in As is 
dominated by the increase in £+, and so R m increases with 
a for a fixed outer boundary condition. 

However, R m decreases with the increase in f e , which 
also shows that these jets are thermally driven. In Fig. Illt l. 
R m is plotted with f e for a — (solid), a — 0.01 (dotted) 
and a — 0.02 (dashed). The fixed parameters are E- = 
0.001 and Ao = 1.65. It was also found out that there are 
critical f e beyond which no standing shock conditions are 
satisfied, and they are f ec = 0.135 for a = (solid), f ec = 
0.22 for a = 0.05 (dotted) and f ec = 0.645 for a = 0.1 
(dashed). It is obvious that f ec increases with the increase 
of a. 

It is to be noted that R m is the fraction of matter which 
is shock heated and ejected out as bipolar outflows. To prop- 
erly understand the role of shock in driving bipolar outflow, 
we consider accretion solutions with the same outer bound- 
ary condition, or the case presented in Fig. I lib . In Fig. 112b ,. 
we plot various quantities across the shock e.g.,, R m (solid), 
£+ (dotted), x s (dashed), and ^4 S (long dashed), as a func- 
tion of ct, and for f e = 0, i.e., corresponding to the solid plot 
of Fig. [T]J- In Fig. [I2p, we plot R m (solid), £+ (dotted), 
Xs (dashed), and ^4 S (long dashed), as a function of a, and 
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shock, which signifies that the mass outflow rate decreases 
with the increase of f e . 
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Figure 13. E — Ao shock parameter space for (a) a = 0.05, and 
each curve is for / = 1, f e = (solid), / > 1, /e = (dotted), 
/ > 1, / e = 0.05 (dashed), and / > 1, / e = 0.2 (long dashed), 
and (b) / > 1, / e = 0.1, and each curve is for a = 0.0 (solid), 
0.05 (dotted), 0.1 (dashed), and 0.15 (long dashed). 



for / e = 0.1, i.e., corresponding to the dashed plot of Fig. 
1 1 lb . Increasing a in accretion with the same outer bound- 
ary condition, will decrease A+ and therefore decrease x s 
(dashed). Decrease in x s implies that the shock is formed 
closer to the black hole, where the viscous dissipation would 
be more i.e.,£+ will be higher (dotted). The jet velocity is 
small at the jet base, so large £+ means hotter flow and 
stronger driving of the outflow i.e. .higher Rm- Figure [T2b 
shows that, as the shock dissipation f e is increased, at cer- 
tain a, £+ is reduced while the decrease in As is marginal, 
this causes Rm initially to decrease with a, but eventually 
starts to increase as £ increases appreciably. 

In Fig. 113b . we plot the E — Ao shock parameter space 
for flows with / = 1, f e = (solid), / > 1, f e = (dot- 
ted), / > 1, f e = 0.05 (dashed), and / > 1 f e = 0.2 (long 
dashed). The viscosity parameter for this figure is a — 0.05. 
It is to be remembered / = 1 means no massloss and / > 1 
means presence of bipolar outflows. Therefore / = 1, f e = 
(solid), implies the shock parameter space for RH shocks 
with no massloss and no dissipation at the shock, while 
/ > 1, f e = (dotted) implies non dissipative shock but 
massloss is present. Similarly, / > 1, fe = 0.05 (dashed) 
and / > 1, fe = 0.2 (long dashed) shock parameter space 
for dissipative shocks of various strength and in presence of 
massloss. The parameter space for standing shock shrinks 
with the increase of f e . In Fig. 113b . we plot the shock pa- 
rameter space for fixed f e = 0.1 and various values of a = 
(solid), a = 0.05 (dotted), a = 0.1 (dashed) and a = 0.15 
(long dashed). We know from Fig. [6l that with the increase 
of a the shock parameter space shrinks, however, in pres- 
ence of massloss, flows with low E seems to show no stand- 
ing shock. In Fig. 113b . we see that although the parameter 
space for shock shrinks but low E flow again exhibit standing 



5 DISCUSSION AND CONCLUSION 

Quasars and micro-quasars may show strong jets, and these 
outflows are correlated with the spectral state of the ob- 
ject. Therefore, quantitative estimate of the generation of 
these outflows are required, especially, the relation between 
the mass outflow rates and the viscosity parameter needs 
to be ascertained, this is because the viscosity will deter- 
mine the disc spectral states. Most of these estimates are 
available for inviscid flow, or special viscosity prescription 
or shock condition, but not for the most general shock con- 
dition, i.e., partially dissipative shock. 

In this paper, our main concern has been to estimate 
the thermally driven bipolar outflows from shocked accre- 
tion discs around black holes. Since accretion disc should 
posses some amount of viscosity, a proper understanding 
of viscous accretion disc is required. As outflows are gen- 
erated from the post-shock disc and high energy photons 
should also emerge from the post-shock disc, an investiga- 
tion of estimating the outflow rate and correlating it with 
the dissipation parameters is a must requirement. In order 
to estimate the outflows correctly, a proper understanding 
of the accretion process has to be undertaken, and we pre- 
sented all possible accretion solutions without massloss in 
section 4.1, including the shocked and shock free accretion 
solutions and the dependence of these solutions on viscos- 
ity parameter. While doing so we have compared solutions 
with same inner and outer boundary conditions, and have 
shown that these differences produces a significant difference 
in interpreting the results. 

For solutions with the same inner boundary conditions 
(i.e., Figs. Eli-p), we found two critical viscosity parameters, 
the first one being the onset of shock (q c i), and the other 
being the one above which standing shock disappears (a cu ) 
and generates a shock free global solution which is wholly 
subsonic except close to the horizon, or the ADAF type solu- 
tion. Simultaneous existence of both a c i and a cu will depend 
on whether Ao is small enough to produce a correspond- 
ing O-type inviscid solution. If inviscid solution is already 
shocked then only a cu will exist, while if the inviscid solution 
is I or W type then neither a c \ nor a cu exists. For solutions 
with the same outer boundary condition, i.e., flow starting 
with same E, Xi„j = AK(iinj) at some injection radius Xinj, 
e.g., Fig. 2£i-e, there would be an additional critical viscos- 
ity parameter ai, which would allow a global solution con- 
necting the horizon and the outer boundary Xmj- We have 
also confirmed that fluids in such cases, x 3 would decrease 
with the increase o f a. The decrease of x„ with the increase 
of a is interesting. IChakrabarti fc Titarchukl (|l995l ) for the 
first time argued that the post-shock disc is the elusive 
Compton cloud, which inverse-Comptonizes the pre-shock 
soft photons to produce power law tail. If the shock remains 
strong then we have the canonical hard state and when the 
shock becomes w eak or disappe a r we ha ve the canonical soft 
state. Moreover, iMolteni et all (|l996bl ) showed that if the 
shock oscillates, it does with a frequency ui ~ x~^, where 
P — 1 — i> 3/2. If the shock oscillates then the hard radiation 
from it would oscillate with the same frequency and could 
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explain the mHz to few tens of Hz QPO observed in stellar 
mass black hole candidates. Outburst phase in microquasars 
starts with low frequency QPOs in hard state, but as the 
source moves to intermediate states the QPO frequency in- 
creases to a maximum and then goes down in the declining 
phase, a fact well explained by approaching oscillati ng x 3 
with the increase of viscosity (|Chakrabarti et oE| [2009l . Our 
steady state model also shows that if every other conditions 
at the outer boundary is same, then x s decreases with the 
increase of a, so we expect with the increase of viscosity the 
shock oscillation too will increase, and therefore the shock 
oscillation model of QPO seems to follow observations. 

We have computed the mass outflow rate from both 
non-dissipative as well as dissipative shocks. The mass out- 
flow rate is always of higher entropy than the pre-shock disc, 
which shows that mass-outflow rate is a natural consequence 
of a shocked accretion disc, a fact readily supported by var- 
ious multi-dimensional simulations. In this connection we 
would like to comment that unlike [Chattopadhvav fc Dasl 
|200^ ; [Das fc Chattopadhvavl (|2008l '). we have corrected the 
jet cross sectional area with the projection effect onto the 
jet streamline. We showed that Rm generally decreases with 
/ e . We also showed that the parameter space is significantly 
modified due to the presence of massloss and dissipation 
at the shock. When dissipative shocks are included we do 
see that the relative mass outflow rate decreases which are 
also to be expected. Increasing dissipation would make the 
shocks weaker, which can be identified with the softer spec- 
tral state, and the acc retion disc in the soft state will give 
weaker or not outflow (|Gallo et. alj|2003l '). 

Comparison of steady shock parameter space Figs. 
I]6l9p . suggests that massloss may trigger an instability. It 
shows that parts of parameter space which produced steady 
shocks in absence of of massloss, do not show steady shocks 
in presence of massloss. This is because the post-shock pres- 
sure gets reduced as it looses mass, and the shock moves 
closer to the black hole (Fig. [7]) to regain the momentum 
balance. However, this is not always possible, and the shock 
may oscillate or get disrupted altogether, therefore, there 
should be a massloss induced instability as well. We have 
also shown that the main driver for the bipolar outflow is the 
energy gained through shock. Interestingly, the mass outflow 
rate generally increase with the increasing viscosity param- 
eter. Since the the shock location also decreases with the 
increasing a for identical outer boundary condition, there is 
a possibility that the QPOs and mass outflow rates be cor- 
related — an issue we will pursue in fully time dependent 
studies. 

It is interesting to note that the dissipation parameter 
used in this work has been assumed constant, however this 
would actually depend on accretion rates and the size of the 
post-shock region. As the accretion rate changes the total 
photons generated would change, and similarly as x s changes 
the fraction of photons intercepted by post-shock disc as 
well as its optical depth would change, this would render / e 
variable. 

It has been observationally established that steady jets 
are observed from black hole candi dates, when the s pec- 
trum of the disc is in the hard state (|Gallo et. alll2003l V In 
our model, the presence of strong shock is similar to the 
hard state. We have shown that with the increasing a, we 
get stronger jet, however evolution of jet states with spec- 



tral states is a time-dependent phenomenon. Moreover, the 
phenomena of QPO and growth or decay of QPO are time 
dependent phenomenon too. Our conjecture that the QPO 
will be correlated with the jet states can only be vindi- 
cated through a fully time dependent study, which is be- 
yond the scope of this paper. Furthermore, since we concen- 
trated on the effect of viscosity on accretion disc and out- 
flows, therefore, magnetic field and other realistic cooling 
processes have been ignored. If cooling processes are con- 
sidered then a direct comparison with the observation will 
be possible. As has been noted that a little bit of magnetic 
field will have an important effect on the dynamics of the 
flow (|Progall2005l ). however transonicity criteria will be im- 
portant too. Because a magnetized flow may posses fast, 
slow or Alfvenic waves, and the number of sonic points may 
increase jTakahashi et al]\20Q($ ). If magnetized flows admit 
shocks, then the shock produced, may be even more robust. 
And this hydrodynamic model might well act as the sim- 
pler version of the magneto fluid model. We are working on 
dynamics of magnetized flows and would be reported else- 
where. 

The concrete conclusions we draw from this paper is 
the following. Viscosity is important and affects the accre- 
tion solutions both quantitatively and qualitatively. Shock 
in accretion can be obtained for fairly high viscosity parame- 
ter. Shocks naturally produces outflows, and for fixed outer 
boundary conditions of the disc, shock location decreases 
but mass outflow rate generally increases. This augurs well 
for the model as this is exactly observed in hard to inter- 
mediate hard spectral transitions. However, in presence of 
dissipative shocks the mass outflow rate decreases. Over all 
we see that Rm may vary between few % to more than 10 
%, although in presence of dissipative shocks, the estimate 
of Rm is about few %. We have also computed the shock pa- 
rameter space for accreting flows without mass loss and dis- 
sipation, with massloss but no dissipation and with massloss 
and dissipation and have shown that the parameter space for 
steady shocks shrinks with increasing dissipation. 
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